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1 Introduction 



The small angle Bhabha scattering (SABS) process is used to measure the luminosity of 
electron-positron colliders. At LEPl an experimental accuracy on the luminosity of baja < 0.1% 
has been reached [1]. However, to obtain the total accuracy, a systematic theoretical error must 
also be added. The accurate determination the SABS cross-section, therefore, directly affects some 
physical values measured at LEPl experiments [2,3]. That is why in recent years a considerable 
attention has been devoted to the Bhabha process [3-11]. The reached accuracy is, however, still 
inadequate. According to these evaluations the theoretical estimates are still incomplete. 

The theoretical calculation of SABS cross-section at LEPl has to cope with two somewhat 
different problems. The first one is the description of an experimental restrictions used for event 
selection in terms of final particles phase space. The second concludes in the writing of matrix 
element squared with required accuracy. There are two approaches for theoretical investigation 
of SABS at LEPl: the approach based on Monte Carlo calculation [3-5,7] and semi-analytical 
one[6,9-ll]. 

The advantage of Monte Carlo method is the possibihty to model different types of detectors 
and event selection [3]. But at this approach one can not use in a strightforward way the exact 
matrix element squared based on essential Feynman diagrams because of infrared divergence. 
Therefore, some additinal procedures (YFS factor exponentiation [12], utihzation of the electron 
structure functions [13]) apply to get rid this problem and to take into account leading contribution 
in the higher orders . It needs to be carefully at this point because of a possibihty of the double 
counting. Any way, up to now the next-to-leading second order correction remains uncertain, and 
this is transparent defect of Monte Carlo approach. 

The advantage of analytical method is the possibility to use the exact matrix element squared. 
The infrared problem in the frame of this approach is solved by usual way taking into account 
virtual, real soft and hard photon emission as well as pair production in every order of perturbation 
theory. Any questions with double counting do not arise at analytical calculations. The defect of 
this method is its low mobility relative the change of an experimental conditions for event selection. 
Nevertheless, the analytical calculations have a great importance because allow to check numerous 
Monte Carlo calculations for different "ideal" detectors. 
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Up to now analytical formulae for SABS cross-section at LEPl are published for the case of 
inclusive event selection (lES) when circular symmetrical detectors record only final electron and 
positron energies [10,11]. These define the first and second order corrections to Born cross-section 
with leading (of the order (aL)" ) and next-to-leading (of the order q;"L"~^ ) accuracy as well 
as third order one with leading only. Just these contributions will have to be computed to reach 
required per mille accuracy for SABS cross-section at LEPl. Note that such accuracy selects only 
coUinear (like two-jets final-state configuration) and semicoUinear (like three-jets one) kinematics. 

In this paper I hst full analytical calculation for lES with wide-narrow angular acceptance. The 
first and second order corrections are derived with next-to-leading accuracy starting from Feynman 
diagrams for two-loops elastic electron-positron scattering, one-loop single photon emission, two 
photon emission and pair production. The third order one is obtained with leading accuracy by 
the help of the electron structure function method. The results for leading second and third order 
corrections in the case of CES are given too. 

The contents of this paper can be outlined as follows. In Section 2 the "observable" cross- 
section Uexp is introdcued with cuts on angles and energies taken into account, and the first order 
correction is obtained. In Section 3 the second order corrections are investigated. These include 
the contributions of the processes of pair (real and virtual) production considered in Subsection 
3.1 and two photons (as well real and virtual) emission. In Subsection 3.2 the correction due 
to one-side two photon emission is considered and in Subsection 3.3 - due to opposite-side one. 
The expression for the second order photonic correction is given in leading approximation only, 
while the next-to-leading conribution to it is written in Appendix A for both symmetrical and 
wide-narrow detectors. The latter does not contain auxiliary infrared parameter. In Section 4 the 
full leading third order correction is derived using the expansion of electron structure functions. 
In Section 5 the numerical results suitable for lES arc presented. The correspondence of obtained 
results with another semi-analytical ones is dicussed in Conclusion. In Appendix B some relations 
are given which have been used in the process of analytical calculations and which will be very 
useful for numerical ones. 
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2 First order correcion 



Let us introduce dimentionless quantity 

1 



where Qf = e^9\ (e is the beam energy and 9i is the minimal angle of the wide detector). The 
" experimetally" measurable cross section Uexp is defined as follows 

a^x, = j dx,dx,ed q, d 6,6, dx^dx^d^q^d^q,^ ' 

where X is undetected final particles, Xi (a;,), qi, (5^) are the energy fraction and the transverse 
component of the momentum of the electron (positron) in the final state. Functions 6? do take 
into account angular cuts while function 6 - cutoff on invariant mass of detected electron and 
positron: 

6^ = ^(^3 - e.)e{e. - Oi) , 6^ = 9(9^ - 9+)9{9+ - ^2) , 6 = 9{xiX2 - x,) , 

e_^\Al^ 0+=''-^. (3) 

xie X2e 
In the case of symmetrical angular acceptance 

^3 

6*2 = 6*1 , 6*3 = 6*4 , p = — > 1 , 



but for wide-narrow one 



9s>9^>92>9,, A = ^ > 1 

PI 



Fof numerical calculation ones usually take 

9i = 0.024 , 93 = 0.058 , ^2 = 0.024 + , ^4 = 0.058 - . 

8 8 

The first order correction Ei includes the contributions of virtual and real soft and hard photon 
emission processes 

Ei = Ev+5 + S^ + Sj^ . (4) 

The contribution due to virtual and real soft photon (with the energy less than Ae, A <^ 1 ) may 
be written as follows ( in this case Xi = X2 = 1, qi + q2 = ^) 

Pi 

E^+5 = 2^/^[2(L-l)lnA + ^L-2], L = In , (5) 



where z — {^2)"^ /Q\ and m is electron mass. 

The second term in r.h.s. of Eq.(4) represents the correction due to hard photon emission by 
the electron. In this case 

X = 7(1 - xi, A;^) , X2 = 1 , + qi + = Q , < Xi < 1 - A . (6) 

It can be derived by integration of the differential cross section of single photon emission over the 
region 



->]_2 

pI< z < pI , x"^ < zi^^ < x^pI , -1 < cosip < 1 , (7) 

Vl 



where (f is the angle between vectors and g^, in the same way as it has been done in [10] 
for symmetrical angular acceptance. But at this passage I would like to indicate the principle 
moments of method used largely to obtain the results of the Section 3 and based on the separate 
calculation of the contributions due to coUinear kinematics and semi-coUinear one [14]. 

In coUinear kinematics emitted photon moves inside the cone within polar angles ^-^ < ^0 ^ 1 
centred along electron momentum direction (initial: k\\pi or final: In semicollinear region 

photon moves outside this cones. Because such distinction no longer has physical meaning, the 
dependence on auxiliary parameter ^0 disappeares in total contribution. This is valid for lES as 
well as for CES. 

Inside coUinear kinematics it needs to keep electron mass in differential cross section 



da — 



1 + 2m^ / I 

1 



Si^i g^ \Si 



dV, 



d^qid^q2d^k , i, n (q\ 

dr = 0^ \pi+p2-k- qi- q2) , (8) 

where q = pi — k — qi, Si = 2{kqi), ti = 2{kpi), s = (2pip2) and Pi(p2) is 4-momentum of initial 
electron (positron). If photon moves inside initial electron cone 

si = x{l - x)eH'i, ti = -m2(l - x)(l + 77), g2 ^ -xh^e^ = -€^$1 , 

dr^—e^7r^x(l-x)dxdr)d9^_, < 77 = ^ < ^ , (9) 



and one can derive after integration relative r] 



Pi 



l + x\ e'e' 2x 
m 



1 — X m? 1 — X 



Q{x'pi - z) . (10) 



The r.h.s.of Eq.(lO) corresponds to the contribution of narrow strip with the width 2y^A(l — x) 
centred around hne z = zi in {z,zi) plane, where A = Oq/Oi. Really, the condition 9^ < 9q for 
initial electron cone may be formulated as follows 

^2 



^/z-^/z^\< A(l 



X) 



-1 < COSip < — 1 + 



"iJZxZ 



fill 



If photon moves inside final electron cone 

^ l^m\\ + 0,^1 = -(1 - x)^e''_ , = -e^el = -e'Ol , 

X 



dT = — 6 TT x(l - x)dxdC^r- , < C < „ 
s x^ 



[12) 



and the integration relative C leads to 



2a^ ''r dz ^ ^ 



1191 



P2 



dx 



i + x\ ehV 
1 — — 

1 — X 



2x 
1-x 



(13) 



The r.h.s. of Eq.(13) corresponds to the contribution of the strip with the width 2-^x^(1 — x)X 
around line Zi = x'^z in plane {zi, z). Really, the condition 0^ < Oq for final electron cone may be 
formulated as | f |< where f — k/uj — / e^, and the last reads as 



r- n ^^ 1 ^ \^x\l - xf - - x^zf 
- x^Jz \< x{l - x)\ -1 < cos(/7 < -IH '- — -££z . 



(14) 



2x-sjzz\ 

Having contributions due to coUinear regions now it needs to find the contribution due to 
semicoUinear ones. Supposing m = in r.h.s. of Eq.(8) the differential cross section suitable for 
this case may be written as follows 

ot'd'{)dzdz\{\ + x^) r 1 



da — 



'kQ\z{z\ — xz) 



X 



dx . 



(15) 



.^1 + 2; + 2^J Z\ZCOSLp Z\ + X^Z + 2x^JZ\Z(X)Si^. 

When integrating the first term into the brackets in r.h.s. of Eq.(15) one must use the restriction 
^-y > ^0 or 

I — |> (1 — x)A , — 1 < cos(/7 < 1 ; 

\m-xf-{^^~^zf 



y/zl — \fz |< (1 — x)X , 1 > cos^p > — 1 + 



2.yJZZl 

while for the integration the second one - the restriction | r | > 6'o or 



(16) 



^/Zl — X\fz \> x{l — X)\ , — 1 < COSip < 1 



Jz[- xsfz\<x{\- x)\ , 1 > cos</7 > -IH 

2xJzz{ 



The integration (15) over the region (16) leads to 



(17) 



2a^ f dz f 1 + 



lor f dz r L- 



dx 



X 



ln- + L2)^r+^3^3 



(18) 



Analogous, the integration of r.h.s. of Eq.(15) over the region (17) gives 



2a'^ r dz f I + x^ z 
(^b = ^^ — / dx[ln^^ + Li 



Qi J z 



p\ 



1 — X V x^A^ 



(19) 



The values Lj which enter into Eqs.(18) and (19) are defined as follows 



Li = ln 



x\z-l){pl-z) 



{x - z){xpl - z) 



, La = In 



[z - X )(x P3 - Z 



x'^{x — z){xpl — z) 



, Ls = In 



[z — x'^){xp\ — z) 



{x - z){x'^pI - z) 



Beside this the following notations for 9- functions are used 



= e{x'pi - z) , er = 1 - or = o{z - x'pi) . 

Thus, the may be represented as the sum of (10), (13), (18) and (19) divided by factor 
Ana'^/Ql or 

^""^l^/S / Y^[i^ + Ot^)iL-i) + K{x,z;ps,l)]dx, (20) 



j(a;) 



Pi 



K{x,z;p3,l) 



(i-xy 



1 + x^ 

Further I will use the short notations for ^-functions, namely 



ef^ = 9{x'p^ -z),9,^ 9{pI - z) , 9i^^^ = 1 - 9f , 9^ = 1-9,. 

One may easy to see that for wide-narrow detectors can be derived from for symmetrical 
ones (see [10]) by the change z-integrations limits 



pI 



dz 



(21) 



Pi 



and the substitution p^ instead of p under integral sign. 

The third term in r.h.s. of Eq.(4) describes photon emission by the positron. It may be derived 
by full analogy with E^ except restrictions on variables z and zi, namely 



1 < Z < p^ , X p2 < Zi < X p^ 



„2 



(22) 



The contribution of the coUinear kinematics {k\\p2 and k\\q2) to single hard photon emission cross 
section corresponds to the integration over the regions inside strips with width 2^/z{l — x)X and 
2-^x^(1 — x)X, respectively. It may be written as follows 



^fe||P2,fe||?2 



2a^ I dz 

of J 

1 X, 



l-A 



X 



2x 



X 



In^ 



2x 



X 



^42 



where 



A — d^^ 



2 ) 



A42 + ^4 — ^2 ■ 



(23) 



(24) 



The contribution of semi-coUinear kinematics may be derived by integration (15), taking into 
account the restrictions (16), (17) and (22). The latters correspond to regions outside narrow 
strips near zi = z and zi = x'^z, respectively. The result is 



2a^ jdz V^l+x^ 



Za-^ r dz /■ i + 



dx 



X 



In ^(A42 + Ag) + L2Ag + (Li - 2 lnx)A42+ 



Lsief^ - ^?^) + L4(^4 - O2) 



(25) 



where 



Li = ln 



Lq = In 



{z- pl)ipl- z)x^ 



{xpl - z){xpl - z) 
z — x^pl){xpl — z 



In 



{z -x^pl){xyl- z) 



La — In 



{x^pI - z){xpl - z) 
The S// is the sum of (23) and (25) divided by 4.7^0^ /Ql 
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2 //n2. 



x'^{xp\ — z){xp2 — z) 

z- pI){xpI - z 



{pi- z){xpl- z) 



(26) 



-'H 



a 
2^ 



J dz 1 - 

J ^ J T 



X 



dx 



X 



(L - 1)(A42 + Ag) + K(x, z; p^, P2) 



(27) 



K = iL^(A42 + Ag) + A42L1 + AgL2 + {ef - 9^^^)Ls + {9, - 62)1, 



1 + x^ 



As one can see the auxiliary parameter ^0 disappears in expressions for and T^h, and large 
logarithm acquires the right appearence. Thus, the separate investigation of contributions due 
to coUinear and semi-coUinear kinematics simplifies the calculations and gives also the dipper 
understanding of underlying physics. The experience of this approach is very important for the 
study of CES when it needs to describe events which belong to electron cluster (or positron one) 
in a different way as compared with events do not. 



The different parts in r.h.s. of Eq.(4) depend on auxiliary infrared paramerter A but the sum 
does not. It has the following form: 



P3 1 

Si = 1^1/ J [-A42 + l{iL- l)Pi(a;)(A42 + AiJ) + ^i?)^^' 



1 



p1 



where 



Pi{x)^\^e{l-x-A) + {2lnA + ^)6{l-x) , A^O. 



f28) 



In order to make the elimination of A -dependence more transparent one can use the following 
relations: 

1 P^{x)dx = - y YT^^^ ' / Pi{x)0f^dx = e'^^"^ I YZ^dx , 



1 VZ/P4^ r, \/5/P2^ „ 

r 1 + 
1 -, 



[ Pi{x)A^£dx = 9^9f^^ [ -^dx - 929i"^^ [ -r^dx , (29) 



where A^J = A42 - aIJ . 

The r.h.s. of Eq.(28) is the full first order QED correction to born SABS cross section at LEPl 
for lES with switched off vacuum polarization. The latter can be taken into account by insertion 
the quantity [1 — Il{—zQ'\)\~'^ under sign of z-integration (for 11 see [3] and references therein). 

3 Second order correction 

The second order corection contains the contributions due to double photons (real and vrtual) 
emission and pair production. As in symmetrical case one needs to distinguish between the 
situations when additional photons attach only one fcrmion line (one-side emission) and two 
fermion lines (opposite-side emission) in corresponding Feynman's diagrams. 



3.1 The contribution of pair production 

Consider at first the contribution of the process of electron-positron pair production E^"'^ to 
the second order correction: 

9 



In order to get rid of the writing some formulae which have the same structure for both symmetrical 
and wide-narrow angular acceptance I will often send the reader to work [11] in which the details 
of computation are given for symmetrical case. 

The experience of Section 2 allows to write the expression for Yf^^~ when created electron- 
positron pair press to electron momentum direction, using the result of [11] for YT^^ suitable for 
wide-wide angular acceptance. It needs only to change z-integration limits; (p^, 1) — > {pi, pi) and 
substitute ps instead of p everywhere under integral sign. The result may be written as follows: 



a 

4^ 



Pi 



dz 



L L l + -ln(l-Xc) - 



J 3 3^ 



X 



__ln(l-a;,) + -ln2(l-a;J+ j -^9^ - - x)j + 



+ 



1 

/ 



LR{x){l + ef^) + erC^{x, z- ps) + C^ix) + d^ix, z- Pa) 



3(0;) 



dx 



(31) 



R{x) = (1 + x)(lnx 



X 



6x 



(4 7x 4x^) , 



^ . , 113 142 2 2 

Ci(x, z\ Pa) = — — + —X - -x 



3x 



-3Lj2(l - x) 



C2{X) 



+ (Sx^ + 3a: - 9 



3' 

^(l + a;)ln(l-a;)+ 
2{5x^ - 6) 



2(1 + x2) 



X) 



2 In 



.2^2 



X'P3 - Z 



Xp3-Z 



X 



X 



.)\nx + ^ ^ In^ X + R{x) In -^^'^^ ^)' 



1 -x 



Pi 



122 133 4 2 2 4, , , 2(1 + x2)^ . . 



+^(-8a;2 - 32x - 20 + 3 + 



13 



2(1 + x^) , 

02(2;, ^; P3) = ^71 r In 

3(1 — X) 



X 1 — X 

{z-x'){pl-z){z-l) 



)lnx + 3{l + x)hi^x, R{x) ^ 2R{x) + -{1 + x) , 



+ +R{x) In 



{z-x'){pI-z){z-1) 



x'^pI — z 



(32) 



{Z - X)2(x2p3 - Z) 

The r.h.s. of Eq.(31) docs not contain infrared auxiliary parameter because it includes the con- 
tributions due to real and virtual pair production. The contribution of hard pair takes into account 
both, coUinear and semi-coUinear kinematics, and this ensures the next-to-leading accuracy. 

If created elctron-positron pair is emitted along of the positron momentum direction the cor- 
responding expression requires more modifications. The source of such modifications is the semi- 
coUinear kinematics as we saw in Section 2 for the single photon emission. 
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The strightforward calculation shows that for contribution of the semi-coUinear region 
(I use here notation p± for 3 - momentum of created positron (electron)) one has to write into 
formula (28) of [11] 



(A42 + AlJ)ln^ + A42ln 



A2 



(^4 - O2) In 



{z-pl){pi-z) 



{z-xpl){xpl- z) 



+ Ai?hn 



{z-x^pI){x^pI- z 



{z- pl){xpl- z) 



[z - xpl){pl- z) 



+ {ef - e^^^) In 



x'^{z — xp2){xpl — z) 

z - x^pl){xpl - z 



+ 



{z - xpl){z - x^pI) 



(33) 



instead of expression in curie brackets and change the upper limit of z-integration: p ^ ps . 

For the contribution of semi-coUinear region the correspnding expression is (see Eq.(33) 

in [11]) 



A42(ln^ + ln 



{z-pI){pI-z) 



+ (^4 - O2) In 



Pliz - pI) 



pUz-pI) 



(34) 



and for semi-coUinear region p_||j?L (see Eq.(38) in [11]) 

(z — x'^pI){x'^pI — z 



All^fln^ + ln 



xjx^plpl 



+ 



pI{z - x'pi) 



For the symmetrical wide-wide angular acceptance pa 



pI{z-x'^pI) 

p , P2 = 1 , and 



(35) 



i42 



e{p' - z)9{z - 1) , Ag ^ e{x'p' - z) , er ^ 9{z - x'p') , , 02, ^ , (36) 



and (33), (34), (35) reduce to corresponding expressions derived in [11] . 

The modification of the contributions due to virtual, real soft and hard coUinear pair production 
includes the change of z-integral upper limit : p — > ps and trivial change of ^—functions under 
integral sign, namely: Qix^p^ — z) — > l^l-l ■, 1 ~^ ^42- The sum of all contributions has the 



following form: 

r dz 



Pi 



A42(l + ^ln(l-xJ)-^/^ 



X 



-i.x) 
^42 



+ A42(---3C2- 



40 



1 



-^ln(l-xc) + ^ln^(l-xe))+/ ^ 



+AgCi(x, z- p2) + A42(C2(x) + d2{x, z- P2)) + {Q^ - 9, 

(xVi - z)pI 



X 42vy-3Mi-^)) + 



^ 9 3 

(-) .Wn^2(1 + x2) 



3(1 -x) 



In 



Li?(a:)(A42 + Ag) + 

{x'^pI- z){xpl- z) 



{x'^pI - z){xpl - z) 



+ 



+R{x) In 



a;2p4 - z)pI 



{xpl- z){z - pI) 



{xpl- z){z- pI) 



+ 
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+R{x) In 



(P2 - ^)pI 



,2\2 



d2{x, z- P2) = -i^^ In — + 2R{x) In 



3(1 - x) (^ - xpi) 



P2 



(37) 



By the help of (36) one can verify that r.h.s. of Eq.(36) goes over in corresponding expression for 
symmetrical angular acceptance. 

3.2 The contribution of one-side double photon emission 

In this Section I give the analytical expressions for all contributions into the second order cor- 
rection which appear due to one-side two photon (real and virtual) emission. The master formula 
which does not contain infrared auxiliary parameter A is written only for leading approximation, 
and next-to-leading contribution to it is given in Apendix A. 

As before it needs to differ the radiation along electron and positron momentum directions 

E2 = E^^ + E^^ , E^^ = E(^+^)' + E(^+^)^ + E^^, 



-"71 



(38) 



The contribution of virtual and real soft photon is the same for both the electron and the 
positron emission 



L(21n2 A + 31nA+ -)- 



8' 



3 4c; 

4WA-71nA + 3C3--C2-- 

2 lb. 



(39) 



Virtual and real soft photon correction to single hard photon emission already differs for 
photon moving along the electron momentum direction and the positron one. In the first case 
corresponding contribution may be derived by the help of result for symmetrical detector (see [10], 
formula(50)) using the substitutions (p|, p|) instead of (p^, 1) for z-integration limits and p^ 
instead of p under integral sign. Therefore, 



„ pI , l-A^ „ 

pI ^= 



+(L-l)(l-F^f ) 



27r2 

pI 

1,2 {I-Xf 

+ - In^ a; - —, ^r- 

2 2 1 + ^2 



/s:(a;,^;p3,l)+ 



+ (l + ^f^)(-2 + lnx-21nA) + ^^' 



-Llnx-|- 



12 



+ 2 In A In X — lna;ln(l — x) — In^x — Li2{l — x) — 



x(l — x) + Axlnx 



2(1 + 

In order to obtain the expression for T:(^s+v)h it needs to change in r.h.s. of Eq.(39): 

i) hmits of z-integration: {pi, pi) (pg, 1) , 



(40) 



ii) Kix, z : p3, 1) ^ K{x, z : p,, p^) ; 9^ ^ , df ^ Ag, 1 ^ A 



42 ■ 



(41) 



The contribution of two hard photons emitted along electron momentum directon may be obtained 
in the same way as E^'^+^)^, using the known result for symmetrical detectors (see [10], Eq.(54)), 
namely: 



47r2 



/ L ( d ( d 

' J z"^ J J ''"xi(l — X — a;i)(l — 

p2 Xc A 



tHH 

Xi){l - Xi^ 



(42) 



,2 X^(X(1 — XljPg — Z)^ 



B = 7(3( ^ + In 



X2(z - l)(pi - Z)(Z - X2)(Z - (1 - Xif )\plx{l - Xi) - Z) 



(p|(i — xiY — zY{z — (1 — xi)y{z — x{i — xi)y{pix'^ — z) 



+ 



C = 7/3(L + 21n 
where 



(l-xi)V 
a;(l — X — xij 

x^plil — XiY — z)^ 



(1 - xif{plx{l - xi) - z){pl{l - Xi) - z) 



-2(l-a;i)/3-2x(l-a;i)7, (43) 



7 = l + (l-Xl)^ /3 = x' + (1 - a;i)2, ( = + (1 - a;l)^ 

7^ = a;a;i(l— x— Xi)— a;^(l— x— Xi)^— 2(1— a;i)/5, 5^ = xxi{l—x—Xi)—x\{l—x—Xif'—2x{l—Xi)'-y. 

Unfortunately, it is impossible to give such simple prescription as (41) in order to obtain T^hh 
from Eqs.(42) and (43). In the case of radiation two hard photons along the positron momentum 
direction an additional detailed consideration of semi-coUinear kinematics is required. All essential 
moments of such consideration shown in Section 2, and reader can make all calculations by the 
help of formulae given in Appendix B of ref.[10]. Here I give final result 

, pI , 1-2A 1-x-A 

/ —^J-^ dx dxi — -— 

J Z'^ J J Xi (1 — X — )( 1 



-'HH 



47r2 



a;i)(l - xi^ 



(44) 



Ihh = AA^S + CA?r^^ + 5A42 + {eP - + - + {04 - ^2)^ , 
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a = 7/3 In 



{plx{l - xi) - z){plx^ - z) 



{P2X{1 -Xi)- Z){plx'^ - Z) 



h = 7/3 In 



{pl{l-xi) - z)(pI- z) 



ipUl-x^)-z){pl-z) 



c = 7/5 In 



{plx{l - xi) - z){pI{1 - xif - zf{pl{l - xi) - z) 



{plx{l - xi) - z){pi{l - x,y - zfiplil - Xi) - z) 



A = 7/5(- + ln 



{p\x^ - z){p\x^ - z) 



5 = 7/3(- + ln 



x^{p\x{\ - xi) - z){p\x{\ - Xi) - z) 

^'^{pI- z){pI- z) 



{pl{l - xi) - z){pl{l - Xi) - z) 



+ Cln 



XX I 

(1 -Xi)^Xi 
x{l — X — Xi) 



+ 5b , 



C = jpl^L + ln 



x'iplil - x,r - zfiplil - x,f - zf 



(1 - - xi) - z){plx{l - Xi) - z){pl{l - xi) - z){pl{l - Xi) - z) 

-2(l-xi)(/? + x7) . 



As one can see the separate contributions in r.h.s. of Eq.(38) depend on infrared auxiliary 
parameter A but E'"''"' and E^^ do not. In order to eliminate A-dependence analytically it needs 
to apply a lot efforts. Below I give leading terms and for next-to-leading ones see Appendix A. 



1 



47r2 



E-^ -^I'^L^Jdx [1(1 + ei^^)P.ix) + I '{P.m (f ) 0: 



2 ''3 , 1 
yL l^dz 

1 Xc 



(t) 

^42 



(45) 



(46) 



where 



P2(x) = Pi®A= / ^Pi(t)P,(^) = lim{ 



(21nA+-)2-4C: 



+ 2 



\-\- 3 1 

(2 ln(l — x) — Inx + -) + -(1 + x)\yix — 1 + X 



5(1 -x)+ 

^(l-a;-A)|, 



(47) 



P2{x)dx = . 

The expressions (45) and (46) are not convenient for numerical calculations. The suitable ones 
may be written as follows 



9 pI , Xc pI , Vz/p3 



pI 



^i(%(|) + ^^2(a;) 



"123 



Xc 



pI 



a 



77 



47r2 I 



p2 mi4 Xc 



P\ , \/^/P4 

az 



dx\^ , (48) 
(49) 
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+ 



where 



g{y) = y + y + 2 ln(l - y) , mas = maxi^p^ , x^Pg) , 

= max{l, ^IpI) > ""^12 = max{l, x^pl) . 

The last two formulae can be derived by means the relations given in Appendix B. The integration 
relative a;- variable in Eqs.(45) and (46) may be performed by the help of the following formulae 

XX X 

J P2{y)dy = F^ix) , j PM9{^)dy = F,{x) , J P,{y)dy = -g{x) , x<l, (50) 



F2{x) ^-2x- — + {x + —) In jYZr^ + 41n(l - x) In - + ALi2{x) , (51) 

2 2 22 

Fg{x) = + {2x + x^) lnx+ {xc + y) In j^-^^ + (2xc + ^ - 2x - y ) ln(x - Xc)+ 

1 — X 



+ 4Li2{x) + 4Li2 



1 — Xr 



Xr < X < 1 . 



(52) 



Therefore, the second order leading contribution to SABS cross section at LEPl can be expressed 
through integral relative z- variable only. 

It is useful to note also that for CES the leading contributions in all orders of perturbation 
theory take into account the emission of photons in initial state only. Thus, the corresponding 
correction due to one-side two photon (real and virtual) emission will be read in this case as 
follows: 



^77 CES 



1 



Pi 1 
pi 

r dz 
~ J ^ 



F2(^]-F2{X,) 
P4 



^(Xc) 



(53) 



(54) 



F2[^]-F2{Xc) 
P2 



3.3 Second order correction due to opposite-side photon emission 



In this Section I calculate analytically the expression for 



^7 — ^S+V + ^S+V + + ^H- 



(55) 
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The quantity does not depend on infrared auxiliary parameter A because it contains all con- 
tributions due to virtual, real soft and hard photon emission. 

The first term in r.h.s. of Eq.(55) takes into account only " oposite-side" virtual and real soft 
photon corrections 



2 ''4 , 
^S+V — ~^ / ~2 



L(41n^A + 61nA+-) -6 - MlnA-Sln^A 
4 



(56) 



The contribution of one-loop virtual and real soft photon corrections to hard single photon 
emission may be written as follows 



^s+v 



2 ^4 , 

a f az 



l-A 



27r^J z 



a r dz 



27r2 J z^ 
1 



2(L - l)lnA + -L - 2 



2(L-l)lnA + -L-2 



1 — X 



2 r 



1 + 0^^^){L-1) + K{x,z;ps,l) 



, (57) 



l-A 



1+X^ 
1 — X 



(A42 + Ag)(L-l)+ir( x,z;p4,P2) 



dx. 

(58) 



In order to find the contribution of two opposite-side hard photon emission into it is 

convenient to use the factorization theorem for differential cross-sections of two-jets processes in 
QCD [16]. It reads as: 



^2 l-A , , , 

Xc 



d i— " ^ ^ 2 1 -|- 

— dxi / dx2- -^xi, z, IPs, 1)^X2, z; p4, P2) , 

J Z J J i — X-\ i — Xo 



(59) 



$(x, z, ; p3, 1) = (A31 + Ai?)(L - 1) + %^(A3i + Ai?) + A31L1 + A^^L^ , 



1 + ^2 

(^f^-^S^))L3 + (^3-^i)ln 



{xpl-z){z-l) 



{z-x){pi-z) 
$(a;,2;, ;p4,p2) = (A42 + ^a2){L - 1) + K{x,z]pi,p2) , 

A3i = ^3-^i, /sS^et^-e^r\ e^^e{i-z), ef^^e{x^-z). 



(60) 
(61) 



The A-dependence of separate terms in r.h.s. of Eq.(55) can be eliminated analytically in 
the whole sum. The leading contribution is expressed in terms of electron structure functions as 
follows 

T 00 1 1 

■ J^L'Jdx, J rfa:2Pi(xi)Pi(a:2)(A3i + Ag^))(A42 + A&V (62) 



E^^ = - 

7 47^2 



Xc Xc/Xl 

The next-to-leading contribution to S!^ is given in Appendix A. 
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The form of suitable for numerical counting may be written in terms of functions F2{x) 
and Fg{x) in the same manner as it was done at the end of Subsection 3.2 



7 4^2 



{-/ 



dz 



i{l)F,{x^) + 2(1)( F,( ^ ) - )^^^ 



P3 



pI 



+ 



XcP3P4 



P4 



+ 



XcP2 
1 



dz 



XcPi 

pI 

f dz 
J 



P4 



XcP3P2 



9 



X. 



(63) 



P3 / "VV^^ VP2/V VPs/ '\^fzJ. 
In the r.h.s. of Eq.(63) the figures into brackets are suitable for CES, when only initial state 
radiation it needs to take into account. 



4 Third order correction 

Inside the required accuracy it needs to keep only leading contribution into the third order 
correction. The latter becomes more important than next-to leading one for LEP2 because of 
increase of the energy. In order to evalulate it one can use the iteration up to the third order of 
the master equation for the electron structure function [13] 



D{x, aeff) = D^'^{x, aeff) + D'^{x, a^ff) . 



(64) 



The iterative form of non-singlet component of Eq.(64) reads 

D^'{x^ a.,,) = 5(1 - X) + E ^(^) V,(x)^^ 

1 

Pi{x) ® ■ -^^ ® Pi{x\ = Piixf, Pi{x) ® Pi{x) = J Pl(t)Pl( 



x\ dt 

IJT 



(65) 
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Up to third order singlet component of Eq.(64) looks as follows [13] 



\ ( Oi 



2! V 27r 



3! V 27r 



2Pi ® R(x) - -R(x) 



(66) 



where R(x) is defined by Eq.(31). Effective couphng Ueff in Eqs. (64) - (66) represents integral 
of running QED constant 

r adt 3, aL\~^ 



a 



eff 

27r 



at/Sn) 2 V 37r/ 



The nonsinglet structure function describes the photon emission and pair production without 
taking into account the identity of final fermions, while singlet one is responsible just on identity 
effects. 

Up to third order the electron structure function has the following form 



D{x,L) = 6{1 -x) + ^P,{x) + \{^)\P2{x) + Ip,{x) + R{, 

, R (x) = Pi © R{x) . 



1 /aL\3rl „ _ „ _ 4 „ _ 2 



, -Psix) + P2(x) + -Pi(x) + -R(x) + R^(x] 

For functions Pz{x) and R\x) see [6,13 MS]. 

The factorization form of the differential cross-section [16] leads to 

oo^ 1 1 

= J ^ J dxi J dx2C{xi,L)C{x2,L) , 



Xc/xi 



Cix„L) = / ^D(t)D(:^)Ag , Cix2,L) = J ^D(t)D(|)Ag . 

Xl XI 



The expansion of C(a;i, L) reads 



aL 



C(xi, L) = 5(1 - xOAg^) + _P,(a;0(Air) + A3i)+ 



C2(a;i)(Ag^) + A3i) + / yAg(:72(xi,t) 



dt 



C3(a;i)(A^-^+A3i) + y yAgC3(a;i,t) 



C^ix) = ^P2(x) + ^Pi(x) + ]^R{x), C2{x,t) = Pi{t)Pi[j 
C,{x) = \p,{x) + \P2{X) + lpi(a;) + ^P(a;) + \r\x) , 



(68) 



(69) 



(70) 
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c,ix,t)^p,it)c4j)+C2{t)p,l^^) , 



(71) 



and the same for C{x2,L) with the substitution X2 instead of Xi and A42^'* (A42) instead of 
Air^ (A31) . 

Because of 9 -functions under integral sign one has to distinguish between 



Sl^jAg and 



In the case of CES one has to acount the initial-state radiation only. Therefore instead of (70) 
it needs to write 



CcEs(x,,L)^A'^, 



(^1) 



(72) 



and analogous for C{x2, L). 

The last step is to write third order contribution in r.h.s. of Eq.(69): 

Xc/x 



3 00 1 i 



(73) 



1 



dt 



Zi = (2A42 + AgAsi + Ag) A42)C3(a:) + J ^(AgAsi + A^S^^2)Cs{x,t) , 

X 



Z2 = [(A31 + AriO(A42 + A'g'') + (A31 + Ari^0(A42 + At^)]Pi{x)C2{x,)+ 

+P,{x) I [AgA42 + AgA3i + Ag)Ag + AgAg]^C2(xi,t) . 

When writing expressions for Zi and Z2 it is taken into account that A31A42 = A42. In the case 
of CES the expressions for Zi and Z2 may be written as follows: 

Z, = (AgA3i + aSa,2)Cs{x) , Z2 = (AgAg^) + A^^^^ Ati)P^{x)C2{x^) . (74) 

Using the relations given in Appendix B the r.h.s. of Eq.(73) may be represented in the form 
suitable for numerical calculations as double integral relative z- and x- variables. It may be written 
as follows: 

E^ = E° + E3 + E^ + E2 , (75) 

where upper (down) index shows the number of additional particles (real and virtual) emitted by 
the electron (the positron). The one-side emission contribute to the r.h.s. of Eq.(75) as 



2 
Pi 



-2 y" Fp{x)dx + 2 J Fr{x)dx- 
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where 



2 

P4 , V^/P4 



-Ot"'' j Fpr{x,Xc)dx - j ^L^O^^"^ j Fpr{x,Xc)dx+ 



+ 



J 



1 Xc 

\/z/P2 

2 ^ Fp^{x^Xc)dx^ , 



r(z, 1) = / R{x)dx = — - + z + z"^ + -z^ - [- + 2z + z'^jlaz , 
J 9 9 V o / 

2 

/(^) = -F,{z) . 



(76) 



In the case of CES the corresponding contribution may be derived by insertion of functions 
Fp, F^ and F^ into the r.h.s of Eq.(76) instead of functions Fp, F^ and Fpr, respectively, where 

F;(x) = C3(X), F;{x) = ^P3(^) + \P2{X) + ^Pl(x) , F;(x) = ^it:(x) + \lf{x) . 

The contribution due to opposite-side emission to r.h.s. of Eq.(75) reads 

2 

A ° 

1 V^/P3 

j R\x)dx -e^^"^ j [H{x,Xc) + 2g{^)h{x;^/p:i))dx 



Xc 



Xc 



, X r 



Y-^L'd^t^ J [^H{x,x,) + 2g{^)h{x-^/p,))dx+ 

1 Xc 



+ 



2 

r37)(a:c) 



j'^L'et^ j {H{x,x:) + 2g{^)h{x-V~z/p2))dx+ 

1 a:c 



2 
P4 



+ 



dz 3 



1 

Vz/P3 



Pi{x)g( 



XcpZPA 
1 



+ 



/ 

Xcpl 



dz 3 



XcPa/\/z 



+^P;^ Mx;^))rfx + (p3- P4) 

\X P4 J \X P4 J \ Pa 



+ 



XcPll\fz 



+ g, 

X P2 / \X P2 



h\x]—jjdx + {p2'r^ 1) 
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2 

dz 3 



XcP3P2 
1 



where 



XcP4 



Xcp4/V^ 



g{a;b)=g{a)-g{b), G{a;b) = G{a) - G{b), G{z) = ^f{z) + ^iz) + ^r{z) , 



If zL If If If 



H{x, Xe) = Pi(x)[2/(^) + -^(^) + r(^; 1)] + ^(-)[P2(x) + R{x)] , 

X ^ X X X 



(77) 



h{x;^z/p)^ J jm)Pi{j)- 



1 + ^' r^^^^JV^/p-^-Ki-^) 



Note that substitutions inside square brackets concern both, hmits of x-integration and expressions 
under x-integral sign. 

In the case of CES the r.h.s. of Eq.(77) requires the following modifications: i) coefficient at 
P-i{x) has to be reduced eight times, coefficients at P2{x) and {x) - four times; ii) it needs to 
suppouse /i = and to substitute H''{x,Xc) instead of H{x,Xc), where 



H\x,Xc) = Pi{x) 



k-) + |^7(-) + k-;i) 

2 X 3 X 2 X 



U^mix) + R{x)] . 



The numerical results 



The numerical calculations carried out for the beam energy e = 46.15GeV, and limited angles 
of circular detectors as given after Eq.(3). The Born cross-section 



47ra^ f dz ( zOj' 



(78) 



Pi 



(in symmetrical wide-wide case the limits of integration are 1 and p|) equals 175.922nb for ww 
angular acceptance and 139.971nb for nn and wn ones. Formula (78) takes into account the 
contributions of the scattered diagram as well as the interference of scattered and annihilation 
ones. The contribution of pure annihilation diagram is proportional to df and is negligible even 
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on the born level. Note, that one has to reduce twice the coefficient at under integral sign in 
the r.h.s. of Eq.(78) if he want restrict himself with the contribution of the scattered diagram 
only. When calculating the QED corrections to the cross-section (78) I systematically ignore 
the terms proportional Of, which have the double logarithmic asymptotic behavior [17] and equal 
parametrically to {a\t\) ln^(|t|/s)/(7rs). The last value is about 0.1 per mille as compared with 
unit for LEPl conditions. 

The results of the numerical calculations of QED correction with the switched off vacuum 
polarization are shown in the Tables 1—3 . For comparsion we give also the corresponding 
numbers derived by the help of Monte Carlo program BLUMI based on the YFS exponentiation 
[3]. 

As one can see from the Tablet there is an approximately constant difference on the level 
of 0.3 per mille between our analytical and MC results inside first order correction. Because 
BLUMI compute the first order correction exactly [18] it may be think that this distiguish is 
caused by omitted in the present calculation terms mentioned above. 



first order correction second order correction 





blumi WW 


WW 


nn 


wn 


blumi WW 


WW 


nn 


wn 


0.1 


166.046 


166.008 


130.813 


134.504 


166.892 


166.958 


131.674 


134.808 


0.3 


164.740 


164.702 


129.797 


133.416 


165.374 


165.447 


130.524 


133.583 


0.5 


162.241 


162.203 


128.001 


131.428 


162.530 


162.574 


128.474 


131.127 


0.7 


155.431 


155.390 


122.922 


125.809 


155.668 


155.597 


123.206 


125.225 


0.9 


134.390 


134.334 


106.478 


107.945 


137.342 


137.153 


108.820 


109.667 



Tablel. The SABS cross-section (in nb) with first and second order photonic correction 

In the Table2 I give the absolute values of the second order correction to SABS cross-section 
taking into account both leading and next-to-leading contributions. The correction due to pair 
production is small in accordance with the results of the work [6]. The second order photonic 
correction is represented as a sum of leading contribution and next-to-leading one. As one can see 
the next-to-leading part is not negligible . 
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pair production two photon emission 





WW 


nn 


wn 


WW 


nn 


wn 


0.1 


0.007 


- 0.004 


0.015 


0.742+0.208 


0.679+0.182 


0.249+0.091 


0.3 


- 0.033 


- 0.033 


- 0.020 


0.546+0.199 


0.556+0.171 


0.069+0.098 


0.5 


- 0.058 


- 0.050 


- 0.041 


0.140+0.231 


0.291+0.182 


- 0.314+0.134 


0.7 


- 0.090 


- 0.074 


- 0.069 


- 0.027+0.234 


0.117+0.187 


- 0.571+0.170 


0.9 


- 0.142 


-0.115 


- 0.115 


2.961-0.142 


2.458-0.116 


1.822-0.090 



Table2. The second order absolute correction to SABS cross-section (in nb) 



In the Tables the absolute value of the leading third order correction and SABS cross-section 
with all corrections obtained in this work are shown. The third order one takes into account three 
photon emission and pair production accompanied by single photon radiation. At large values of 
Xc this correction is comparable with second order next-to-leading one. This effect will increase 
in the conditions of LEP2. 



third order correction SABS cross-section at LEPl 





WW 


nn 


wn 


WW 


nn 


wn 


0.1 


- 0.055 


- 0.047 


- 0.006 


166.910 


131.623 


134.817 


0.3 


- 0.065 


- 0.053 


- 0.018 


165.349 


10.438 


133.545 


0.5 


- 0.036 


- 0.040 


0.004 


162.472 


128.384 


131.090 


0.7 


0.089 


0.058 


0.124 


155.596 


123.190 


125.310 


0.9 


0.291 


0.220 


0.331 


137.307 


108.927 


109.893 



Tables. Leading third order correction and SABS cross-section as obtained in this work 



As concerns the second order correction it needs to have the analytical formulae based on expo- 
nentiated form of the electron structure function in order to be consequent in the comparison with 
the BLUMI results. On the other hand, the comparison of given here the second order photonic 
correction, which includes the leading and next-to-leading contributions, with the corresponding 
numbers for non-exponentiated BLUMI version [3] was done recently in [22], and the agreement 
is very impressive. 

6 Conclusion 

In this paper analytical calculation of QED correction to SABS cross section at LEPl are 
given for the case of inclusive event selection and wide-narrow angular acceptance. These include 
leading and next-to-leading contributions in first and second orders of perturbation theory and 
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leading one in the third order. The leading contributions in the case of calorimeter event selection 
are obtained too for any form of final electron and positron clusters. Results are represented in 
the form of manifold integrals with definite limits, and functions under integral sign have not any 
physical singularities. No problem arises with infrared divergence and double counting. 

The selection of essential Feynman diagrams, utilization of natural for this problem Sudakov's 
variables, impact factor representation of differential cross section due to t-channel photon ex- 
change as well as electron structure function method and investigation of underlying kinematics 
were very useful along of the whole this work. It needs to emphasize separately the simplifica- 
tions connected with impact factor representation which allows to represent the differential cross 
sections of two-jets processes in QED by factorized form. The latter allows to use cut-off 9 func- 
tions for the final electron and positron independently on the level of the differential cross-section. 
The calculation does not require to go to c.m.s. of underlying subprocess (as in [6]) and escapes 
corresponding complications. 

At this point I want to comment the analytical calculation of leading contribution due to 
photon emission and pair production carried out in [6]. Authors of these articles used as the 
master formula for description QED corrections to the SABS cross-section due to initial-state 
radiation the representation vahd for cross sections of Drell-Yan process [19], electron-positron 
annihilation into muons (or hadrons) [20] and large angle Bhabha scattering [21]. But inside this 
set the SABS process has a very particular feature caused by the existence of two different scales. 
The first one is the momentum transfer squared t, and just this scale defines the value of the cross- 
section. The second scale is full c.m.s. energy squared s = 4e^, and the quantity 9^ ~ \t\/s « 1 
has status of a small correction. 

The t-scale physics is very simple and defined by peripheral interaction of the electron and 
the positron due to one photon exchange, provided momentum transfer is pure perpendicular : 
t — — g^. The s-scale physics is more complicated. On the born level it exhibits as contribution of 
an annihilation diagram and beside this permits the energy and longitudinal momentum exchange 
for the contribution of scattering diagram. The first order QED correction for s-scale cross-section 
includes the contributions of box diagrams, large angle photon emission and up-down interference 
because both, the eikonal representation for the scattering amplitude and the factorization form 
of the differential cross-section, breaks down. In the second order large angle pair production and 

24 



appear. 



The structure function used in [6] controls i-scale cross-section only and has not any relation 
to s-scale one because physics of different scales evolute by its own laws. 

On the other hand, only scattered diagram contributes in born cross-section used in [6]. But 
everytime when somebody neglects annihilation diagram as compared with scattering one he must 
automatically neglect 6*^ as compared with unit everywhere including the born cross-section (see 
comments to Eq.(78)) and experimental cuts in order to be consequent. Taking into account 
these arguments the master formula in [6] must be necessary simplified by eliminating terms 
proportional ^ ~ \t\/s « 1 and in the numerator of Eq.(5) and in the cutoff restrictions. After 
this it becames adequate to one obtained in [10] and used in this work. 

Numerical evaluations shows good agreement with Monte Carlo calculations inside first order 
but the achievement of an agreement for high order corrections will require an additional efforts, 
connected with writing the version, based on the exponentiated form of the electron structure 
function for present analytical calculation. 
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Pietrzyk for fruitful discussions and critical remarks as well as A. Arbuzov and G. Gach for 
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Let us begin with the consideration of the next-to- leading second order A-independent contribution 
due to one-side two photons emission. At first I will give analytical expression for symmetrical 
case, because it was not published up to now in relevant form. (1 do not introduce special notation 
for next-to-leading contribution to E keeping in mind that only such kind terms are considered 
along this Appendix) 
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(A.l) 



1 



1 l-x 




dx <A + 



dxi — 
Ixi 
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For wide-narrow angular acceptance it needs to consider only the case of the positron emission 
E^^, because the corresponding expression for the electron emission E^^ is just eq.(Al) with (p|, p|) 
as the limits of ^-integration and p3 instead of p under the integral sign. 
The analytical expression for E-y^ has the following form: 



^-^^"4^7 1 z^ 
1 

= yA42 + j dxj A42 4(4 + 3x) + 6(x + 3) ln(l - x) + (x - 1+ 



(A.2) 



26 



+4 



1-x 



A 



(x) 
42 



(1 -x)(3 + lnx) + 2(a; + 3)ln(l -x) +4 4^^1nx 



+(1 + xf) ln(l - x) + 2 ii±ii! (e.e^f In 



1 — x 



P4- Vz 



1 — X 
- ^2^?^ In 



\/z-Xp2 



P2- 



1-X 



+ 



9 

B+ / dxi 





1-x 



1 + x^ 
(1 — 



A1?/i+ + A42/2+ + (C' - CO^l- + (^4 - ^2)/2- ) + 



+ -1- 



1 +X 



X 



1-xi {l-xiY 
+Ai^,-^^)f21n 



(^)A. (1 -xO^xs 



AfsMln 



+ ^3+ 1 + A42(ln 



(1 — XiYxxi 



X2 



1 + X 



+ /5+ ) + (^f - e'^^%_ + (^4 - ^2)^4-+ 



4 , x 

— m 



X2 1 - Xi 



1 — Xi 

Vi/p4-a; 



^(1-^1) 



^42 



2 1±^»/;' 

1 — a; 



l-x 



+ 



[ ^I, 

^ Xi 



1 — X 



1-X ^ 

/ ^"^(^ 



rfXi ( Iq- 

^Xi 



k+ 



+I1 - )+ / ^£irJ|, 

^2 1 — Xi/ J Xi i j 



B = A42(-21n^x + 21n(l -x)ln 



x^i^ - Plf{z - x^pI){x^pI - z){pI - zf 



[z — xplY^xpl — zY 



+ 



+Ag(ln^x + 21n(l -x)ln 



{z-x^pI){x^pI - 2;) 



X'^{z — Xp2){xp\ — Z) 



+AiJ(7-21nxln(l-x) -21n2x-2Li2(l-x) - 
+2(^4 - O2) ln(l - x) In 



+ (3 - 2\nx)K{x, z] Pa, p2)+ 
x{l — x) + 4a;lna;' 



1 + X2 

{xpl- zf{z - plf{z - x'^pI) 



+ 



(pI - zy{x'^p\ - z){z - xplf 

Z — x'^P2){xp\ — Z 



+ 



(xVI - z){xpl - z) 



+4^4^?^ In 



XpA - \fz 



In 



x{x(?i — z) 



P4 

h± = (1 ± c) In 



x'^pl — z 



+ 4029^ In 



\fz -XP2 



P2 



+ 



In 



z — x'^pI 



x{z — XP2) 



{z ~ x^pI){z - xpl) 



l2± = (1±C) 



In 



(z - x{l - xi)pl){z - x{x + xi)pl) 

{Z - Xplf{z - (1 - Xifplf{z - (,T + Xifplf 



{z - x'^pI){z - x{l - xi)pl){z - x{x + xi)pl){z - (1 - xi)piy{z - {x + Xi)pIY 



k± = (l±c)ln 



z — x'^pI 



z-x{l- Xi)pl 
/5± = (l±c)ln 



, /4_|_ = (1 ± c) In 
{z-{l- x.fplY 



z-pI 



z-{l- Xi)pl 



{z - x{l - xi)pl){z - (1 - xi)pl) 



27 



/fi = In 



/7 = ln 



(1 - X^nz - X{1 - X^)pl)\z - (1 - Xi)pi)2 

{z-xplf{z-{x + xifplY 



{z — x'^piy^z — x{x + Xi)p2y{z — (a; + xi)^!)^ 
where X2 = I — x — Xi, and c is the operator of the substitution 



C/(P2) = /(P4) 
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One can verify that in the symmetrical hmit formula ( |A.2| ) coincides with ( A.l ) one. 

For opposite-side emission the next-to-leading contribution to S in the symmetrical case reads 
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We introduce the following reduced notation for 6'-functions: 
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. The quantity K{x^ z\ p, 1) entering into espression for A is the K- factor for single photon 
emission, and the quantity A^(x, 2;;p, 1) may be derived by the help of Eq.(lO) in the following 
way: 
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In the wide-narrow angular acceptance the corresponding formula for may be written as 
follows: 
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It is obvious that in symmetrical hmit formula (A.7) coinsides with (A. 4) one. 
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Appendix B 



Here I give some relations which were used in the process of analytical calculations and at numerical 
computations. For the case of emission along the electron momentum direction they reads 

pI 1 pi v^/P3 

j dz j dx = Jdz J dx , 

p\ 1 \-x P\ \f^/p?, 1-x 

j dz j dx j dxif^""'^ ^ j dze^"^ j dx j dxi. (B.l) 
For the case of the emission along the positron momentum direction: 

pi 1 pi 1 

jdz Jdx [^f ^ - ^] ^ Jdz Jdx\04-e2 + 040^^ + 020^] 

1 Xc 1 Xc 

pi 1 V^/P4 v/^/P2 

^Jdz^^(04-02)jdx+0404^^ j dx + 0202 ^^ j dx^ (B.2) 

1 Xc 3jc 3jc 

pi 1 1-x Pa 1 1-x 

Jdz Jdx J dxi \0'l~'''^ - 0i^~'''^] ^ Jdz Jdx l^{04- 02) + J dxi 

1 Xc 1 Xc 

V^/P4, 1-x V^/p2 1-x 

+^f'^^4 J dx J dxi + 0'2^^02 J dx J rfxij. 



l-yG/p4 Xc l-^/p2 

Some additional relations arise for the case of the opposite-side emission. Let us consider first the 
integration limits restrictions for the product of ^-functions in the symmetrical case: 

t/st/g t/g , yi^g ^ , Ultfi Ui . (.tS-jj 

At first it needs to use the formulae (B.l) and get rid ^-^^^ using the following changes: i) ^^-^^^ — > 
^{xc/xi) ^ ii) the upper limit of X2 integration in the case of 0^^^'' has to be replaced by {y/z/ps) 
and in the case of 0^^^^ by ^/z. 

Thus, there are three regions defined by following curves in {z,Xi) plane: 
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It easy to see that the hmits of integrations may be transformed as follows: 



1 yz/P Vz 



^cP Xc/y/z Xc/xi 

and for / 6i6j^ 6^ the formulae may be obtained from the above ones by putting p — 1. For the 
wide-narrow angular acceptance the prescription is similar: 

The another variants of restrictions in wide-narrow ansular acceptancee may be transformed as 
follows: 
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1 ^ Vz/P4 

^cP4 XcPi/Vz Xc/xi 
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